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Abstract 

The flat inflationary dust universe with matter creation proposed by 
Prigogine and coworkers is generalized and its dynamical properties 
are reexamined. It is shown that the starting point of these models 
depends critically on a dimensionless parameter S, closely related to 
the matter creation rate tp. For £ bigger or smaller than unity flat 
universes can emerge, respectively, either like a Big-Bang FRW sin- 
gularity or as a Minkowski space-time at t = — oo. The case £ = 1 
corresponds to a de Sitter-type solution, a fixed point in the phase di- 
agram of the system, supported by the matter creation process. The 
curvature effects have also been investigated. The inflating de Sitter 
is a universal attractor for all expanding solutions regardless of the 
initial conditions as well as of the curvature parameter. 
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1 Introduction 



Much effort has been spent to understand the effects of the matter 
creation process on the universe evolution [1-7]. In the framework of the 
Friedman- Robertson- Walker (FRW) geometries, a self-consistent phenomeno- 
logical description for matter creation has recently been proposed by Pri- 
gogine and coworkers || |9|. The leitmotiv of this approach is that the matter 
creation process, at the expense of the gravitational field, can happen only as 
an irreversible process constrained by the usual requirements of nonequilib- 
rium thermodynamics. The crucial ingredient, however, is the explicit use of 
a balance equation to the number density of created particles in addition to 
the Einstein field equations(EFE). When properly combined with the second 
law of thermodynamics, such an equation leads naturally to a reinterpreta- 
tion of the stress tensor, corresponding to an additional pressure term, which 
in turn, depends on the matter creation rate. 

An extended manifestly covariant version of such a formulation has also 



appeared in the literature |L0], [TT[| . The nonequivalence between the mat- 
ter creation process and the mechanism of bulk viscosity, which has been 
widely used in the literature as a phenomenological description to the for- 
mer, has been recently clarified |T2|, 13]. In connection to this, we remark 
that the general thermodynamic properties of the formulation by Prigogine 
and coworkers has been more carefully investigated than its dynamic coun- 
terpart, which was also presented in the mentioned papers as a new example 
of nontraditional cosmology. 
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In this article we focus our attention on the cosmological scenario pro- 
posed in Refs.[§, f|. As will be seen, our study will provide some general 
results often useful to bring forward certain subtleties in the cosmological 
solutions which apparently were not perceived at first by the mentioned au- 
thors. To be more specific, in the scenario proposed in the quoted papers, 
the spacetime starts from a Minkowski phase at t — 0, with a particle num- 
ber density n a describing the initial fluctuation (Cf., for instance, ||, pg. 
773). As will be shown here, however, there is no initial fluctuation since the 
Minkowski spacetime starts at t = — oo with a number of particles precisely 
equal to zero, as it should be. Subsequently, these models evolve smoothly 
to a de Sitter phase, in such a way that at t — all solutions describe an 
expanding (inflationary) flat FRW-type universe. Another unnoticed aspect 
is related to the existence of a new large class of solutions starting from a 
big-bang FRW singularity and also approaching de Sitter spacetime for late 
cosmological times. For all values of the parameter of the equation of state 
p = (7 — l)p , the starting point of the flat universes depend critically on a 
dimensionless parameter closely related to the matter creation rate. It can 
also be shown that all expanding solutions converge to a flat inflating de Sit- 
ter solution. In other words, the latter is an attractor independent of initial 
conditions. 

Using phase space portrait techniques, the above analysis is extended 
to curved spacetimes. Due to degeneracy in flat space, the class of models 
proposed by Prigogine et al. is shown to split into two distinct classes. In 
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particular, it is shown that only de Sitter and a space-like singularity (Big 
Crunch) are attracting nodes (stable points) in the phase space of solutions. 
For the first class, the former attractor is independent of the initial conditions, 
that is, all expanding solutions end up like a inflating de Sitter spacetime 
regardless of the values of the curvature parameter. 

This paper is organized as follows: in Section 2, the basic equations de- 
scribing a FRW-type Cosmology with matter creation are presented. In 
Section 3 we discuss the flat case for a matter content satisfying the 7- law 
equation of state. The flat dust model proposed by Prigogine et al. is dis- 
cussed in detail and a second set of exact FRW-type solutions not foreseen 
by those authors are established. Finally, in Section 4 we develop a qual- 
itative analysis in order to investigate the curvature effects. The matter 
creation ansatz of Refs. M, 5] is generalized and its dynamic consequences are 
discussed. 

2 Basic Equations 

Consider now the homogeneous and isotropic FRW line element 

ds 2 = d ti-Ri(t)(-^— + r 2 d6 2 + r 2 sin 2 (6)d(J) 2 ) , (1) 
1 kt 

where R is the scale factor and k — 0, ±1 is the curvature parameter. 
Throughout we use units such that c = 1. 

In that background, the basic dynamic equations for a self-gravitating 



perfect fluid endowed with matter creation reduce to [11, 12: 
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i) Einstein's field equations(EFE): 



R 2 k 

W = 3^ + 3^ , (2) 



x(p + P c ) = -2^-^-^ , (3) 



whereas the energy conservation law, 



p + 3(p + p + p c )# = , (4) 

is the usual consistency condition contained in the independent EFE (|2|) and 
(H). In the above equations, x — , p, p and p c = —aip/3H are, respec- 
tively, the energy density, equilibrium and creation pressures. An overdot 
means time derivative and H = R/R is the Hubble parameter. 

ii) The balance equation for the particle number density: 



n + 3nH = ip 



(5) 



where ip is the particle source (ip > 0) or sink (ip < 0). 

The phenomenological parameter a appearing in the definition of the 



creation pressure was originally introduced in Refs. [10 , |i~Tf . As shown 
there, the Prigogine et al. formulation is recovered in the "adiabatic limit" , 
namely when if) is different from zero but the specific entropy (per particle) 
is constant. From now on, we consider just the "adiabatic case" for which 
the a parameter is given by (see Refs.|TT|, |T2|| ) 
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P + P fR\ 

a = , (6) 

n 

with the creation pressure assuming the following form 

In order to obtain a definite cosmological scenario with matter creation, 
one still needs to provide two additional relations: the equation of state and 
the matter creation rate. In the cosmological domain, the former is usually 
expressed in terms of the "gamma-law" equation of state, 

p=(7-l)p , 0< 7 <2 , (8) 

where 7 is the "adiabatic index" . 

Note that, from @), @, (0) and the evolution equation for the scale 
factor R can be written as 

v 2 2nH ! v 2 2nH ! ' v ' 

which depends only on the matter creation rate if} and the number density 
of particles n. 

Furthermore, using (^) and (0), it is easy to show that n and p satisfy 
the general relation (see Ref. ||15|| ) 



n = n (^y , (10) 

Po 
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so that, using the EFE (fj) one can obtain n(R,H 2 ), with Eq.(^) now de- 
pending only on ip. 

In what follows, we will first discuss the simple and exactly integrable flat 
case (k = 0). This class includes Prigogine et aVs dust model || [|. The 
contributions of curvature terms will be analyze in section 4, but for now 
it suffices to say that flat solutions are not necessarily stable with respect 
to curvature perturbations. In the language of dynamical systems, FRW 
cosmology without matter creation is "structurally stable" with respect to 
curvature perturbations, while models with matter creation are not. 

3 Flat Case: Two Classes of Models 

Let us now discuss the flat case. In terms of the Hubble parameter the 
evolution equation @ for the scale factor reduces to: 



2 2n K ' 

while from (0) and ([H]) the particle number density n assume the form below, 



n = n (^-) 2 - . (12) 

Now, in order to properly generalize the dust universe of Refs.|§, § we 
consider a class of models endowed with the same matter creation rate 



ip = 3(3H 2 . (13) 
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In the above expression (3 is a new phenomenological parameter in the 
notation of Refs.|| ||). Parenthetically, since we are dealing with a spa- 
tially flat spacetime, the above equation implies that ip is proportional to 
the energy density, more precisely, ip = x@P- Of course, these relations are 
not equivalent in curved spacetimes, and so doing we are effectively working 
with two different creation rates (see section 4). 
Inserting Eqs.(|T2|) and ( |1"3"D into ([[1]) one obtains 
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H 

wj 1 



( 7 ^2) , (14) 



and 



H + 3(1 - (^)H 2 = ( 7 = 2) . (15) 
In Eqs.(P~2"D-(|H)D, n is the value of n for H = H Q and H d is given by 



H d = H (-^r 2 (7/2) . (16) 
The case 7 = 2 is the simplest one to be analyzed. However, since it does 



not exhibit the peculiarities present in the class of solutions of (|14|), which 
contains the model of Refs. ||, |J, we do not consider it in what follows. 
Apart from some interpretation problems, Prigogine's et al. model must be 
recovered by taking 7 = 1 in the general solution of (|I~4"D. A straightforward 
integration of this later equation yields 



3(-y_2) 7 

H = H d (l + AR*)^ , (17) 
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where A is a dimensional integration constant. 

Now, in order to introduce in our approach a useful dimensionless param- 
eter, we compute the value of A as a function of H d , H Q and R a , where R a 
is the value of R when H = H Q . From fllTI) one reads 



A = (E - 



1(2-7) 
o 



where 



In terms of E the Hubble parameter takes the form 



H = H d 
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(18) 



(19) 



(20) 



and from fll2|) and (p!7|)-(^9|), the particle number density is given by 



n = n n Tji- 2 



2 

2- 7 



(21) 



with the time-dependent net number of particles, iV = nR 3 , assuming the 
following form: 



N(t) = n Q R 3 ET- 



1-(E-1)A§M 

JX 
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(22) 



For these models, the qualitative behavior at early and late times can 
be easily determined from Eqs.([20|) and (pl[) or (|22|). If E = 1, the Hubble 
parameter is constant, H = H d , with the solutions reducing to a flat de 
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Sitter-type universe regardless of the value of 7. In this way, we may have, for 
instance, a de Sitter universe supported by the adiabatic creation of photons 
(7 = 4/3) or dust (7 = 1). Such a fact had already been observed in Ref.[|]. 
Indeed, H = Ha is a special solution, a singular point in the phase diagram 
of the system {H(p,H), p(p,H)} (see Eqs.(|3~4]) and (|35| ) in section 4). It 
should be noticed that for large values of R, that is, (R/ R )~ 3( - 2 ~' y ^ 2 <C 1, 
the models evolve to this inflating solution. 

How these flat universes emerge depends only on the sign of A or equiva- 
lently, if the dimensionless parameter E is bigger, smaller or equal unity. For 

S > 1, the second term on the rhs of (^0[) is dominant for small values of 

2-7 

R, that is, (R/R ) 1 ^> 1, thereby leading to the usual FRW singularity, 
R ~ t 2/37 . Since H -> 00 in this limit, it follows from (0) and (03) that 
p and n are infinite as t goes to zero. These results characterize a class of 
singular solutions not previously perceived in Refs. |8|, |9| . 

On the other hand, if E < 1 it follows from ( p0|) that there is a minimal 
value of R, namely: 

Rmin = Ro{l ~ S)^T , (23) 

for which H = 0. In addition, Eqs.flU]), (U) and (0), yield p = n = N = 
at R = Rmin, making explicit that such models start, for all values of 7, 
as a Minkowski vacuum. We remark that, from (0) and (|T2]), the creation 
pressure p c scales with H Zl ~ 2 . Thus, according to the above results, in the 
beginning of the universe the creation pressure was either zero or infinite if, 
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respectively, E > 1 or E < 1. Only in the de Sitter case (E = 1) the creation 
pressure has a finite(and constant) value. 

The above results are rather important in what follows. For R = R one 
obtains H = H Q , n = n Q and iV = n R 3 for all values of 7 and S ^ 1. For 
E = 1, (|20|) and ( ^H) reduce to the constant values H = Hd and n = n Q 



as required by the symmetries of the de Sitter spacetime. The latter case is 
different from the de Sitter phase attained in the course of the evolution for 
E ^ 1. In fact, for R > R Q , (|l|) yields 
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n d = n E 2 -t , (24) 



which may also be obtained from (|I~2"D taking H = H^- The important point 



to keep in mind here is that models starting as Minkowski (E < 1) have 
n = when R = R m i n , that is, n Q cannot be the number density at the 
beginning of the universe. 

Therefore, since for 7 = 1 the energy density is p = nM (where M is the 
mass of the created dust particles) we can say that the authors of Refs. || |J 
studied only the case E < 1 , that is, M > ^jf. Note that n and (3 are 
coupled to give the natural mass scale of the models so that for 7 = 1 the 
condition E < 1 (E > 1) assumes the interesting form M > M c (M < M c ), 
where the critical mass is M c = , just the mass of the particles created 
in the de Sitter spacetime. For this steady state scenario, we see from (5) 
and (13) that (3 = ridHd" 1 , and using the present data it is readily obtained 
M c « lGev i.e., the proton mass. Naturally, the above condition may also be 
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translated as a constraint on f3, the free parameter of the models. In terms 
of the mass M, the critical value of this parameter is given by (3 C = ^g^(~]^)i 
where M p i is the Planck mass. In contrast with the approach developed 
here, the mass of the created particles in Refs.|| |J has been estimated by 
assuming that the de Sitter phase is unstable and evolves continuously (up 
to first derivatives) to the FRW radiation phase. 

3.1 Prigogine et aL's Model (S < 1, 7 = 1) 

Let us now analyze with more detail the case 7 = 1. As remarked above, in 
this case xPo = xn Q M = 3H% , so that (|T6| ) and ( [24]) reduce to 



H^-f- , (25) 

and 



n d = ^~ ■ ( 26 ) 

The aforementioned results are presented, respectively, in Eqs.(20) and 
(4) of Refs. 0, Q . Now, in order to relate ( |2^| ) and ( |2lf ) with the characteristic 
time scale of the de Sitter phase we insert ( p5l) into ( p0|) and integrate it with 
7 = 1 and S < 1 to obtain 

J R(t)= J R D [l + S(e 2 ^'-l)] 2/3 , (27) 

which is the same expression for the scale function presented in Refs.|§, f|. 
Therefore, after a characteristic time r c = = \Hj l , the universe reaches 
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the de Sitter phase expanding as 



R(t) ~ R^f^eHat 9 (28) 



where rid and Hd were defined by ( pq) and (|26|) . On the other hand, it follows 



from (^3|) that the minimal value of R for 7 = 1 is given by 

R min = R [l-n 2 / 3 , (29) 

which depends on n Q through E. Hence, if we consider R Q = 1 as in Refs.[|], 
|J, one has -R m i„ < 1, with fl29|) making explicit the constraint S < 1 as 
noticed earlier. Of course, the natural choice is to take not R Q but R m in = 1- 
However, the important step is to determine at what time the scale factor 
assumes its minimal value R m in- From fl27|) we see that R = R m in only if 
t = — 00, which is coherent with our qualitative analysis. Note that for t = 
( P7[) yields R = R a and, (JH]), n = n . In this context we remark that the 
authors of Refs.|| 3 truncated arbitrarily the time coordinate at t — 0, and 
since the arbitrary value of R a has been fixed equal to unity, they obtained 
-R(O) = 1 ; in their words: "the universe emerges without singularity at 
t = 0, with a particle number density n a describing the initial Minkowskian 
fluctuation" (see Ref.||, pg.773). Indeed, for t — 0, in addition to R = R Q 



and n = n , one obtains from (0) or equivalently from ([lq ), (|19D and 



( 20|) that the Hubble parameter itself is H = H Q , making it clear that the 
spacetime is not Minkowski. It thus follows that at t = the solution ([27]) 
describes an expanding FRW-type universe driven by the matter creation 
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process. Indeed, such a result is valid for all values of 7 in the considered 
interval and S ^ 1. As a matter of fact, the extension of the time coordinate 
for t = —00 may also be justified looking for the expression of the total 
number of particles given by (0). By taking 7 = 1 in the later equation and 
replacing the expression of R given in (|27D one obtains 

N(t)=N e (3Mt , (30) 

which is the same expression derived in Refs. || || using a different approach. 
Notice that the Minkowskian limit, N — > 0, is recovered only if t —>■ — 00 in 
accordance with our previous comments, while for t = or equivalently 
R = R Q one has N = N a , as it should be. Despite the fact that this "initial 
fluctuation" n Q does not exist, it will be shown that there is a structural 
instability related with these types of models, which may be discussed using 
the dynamical systems technique (see section 4). 

Summarizing, this class of spacetimes (£ < 1), which includes the dust 
case, emerges at t = —00 as a true Minkowski spacetime, that is, p = n = 
H = and R = R m in, with a general expression given by (|23|). There is no 
initial fluctuation. Due to the matter creation process, the universe evolves 
smoothly from Minkowski to a de Sitter spacetime. This is an interesting 
example of spacetimes which are unbounded in time, that is, their evolution 
ranges the time interval (—00, 00) and are also free of physical singularities. 
Naturally, the characteristic quantities of the de Sitter phase may, for certain 
special values of 7, be independent of n ; a particular value of n arbitrarily 
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fixed at t = 0. This is not remarkable and happens just for 7 = 1 as a 
consequence of the relation p = nM (see Eqs.(p5|), ( p6|) and the definition of 



Tr . 



3.2 A New Model (E > 1, 7 = 1) 

Having in mind that the case £ = 1 is trivial we now analyze the last possi- 
bility, that is, £ > 1 or xfi 2 M < 3n Q . In this expected, Eq. (|20|) is 
not modified, but on the other hand there is no R m in 1 so that R — > and 
H — ► 00, characterizing an initial FRW-type singularity. Taking into account 
such remarks in the integration process of (|20|), it is easily found, for 7=1, 

R(t) = i? [(£-l)(e§^'-l)] 2/3 . (31) 

This expression is consistent with our earlier qualitative analysis. For 
times t > r c = \H^ 1 , the model approaches the de Sitter regime. If t <ti r c 
then ([JI]) yields 



R(t)~R [(i:-l)(l + ^H d t + ----l)} 2/3 , (32) 

so that 



R(t) ~ R [^H d (E - l)i] 2 / 3 . (33) 

This was expected, since for large values of £ = H Q / Hd, the equation 
above approaches R ~ R (^H t) 2 ^ 3 , which is the standard form of the FRW 
dust model [TJ . 
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4 General Case: Qualitative Analysis 



The flat case have already showed striking richness, and a natural extension 
to curved spacetimes is already more than justified. However, since is rather 
difficult to analytically derive the complete set of solutions, we will employ a 
different approach provided by the qualitative analysis of dynamical systems 
H], [17|]. The basic idea is to reduce the field equations to a bidimensional 



autonomous system and perform the qualitative analysis of all solutions. As 
we know, the first step is set up a convenient set of dynamic variables, which 
we choose as being the energy density p and the Hubble parameter H. Using 
the definition H = R/R and the constraint (|2]) we rewrite (|9]) as 

^ = i<3^- 1)w+ 5 to - 3ff2) ■ (34) 

The second dynamical equation is just the energy conservation law(||). 
Inserting (ffi) and (||) it takes the following form 

^ 3 ^- 1) " ff ■ (35) 

Since the ratio \[>/n is a function of the energy density or Hubble parameter 
(or some combination of them) our dynamical system is fully defined by 
Eqs.(^4[) and (0). In what follows we consider only a positive-definite 7. 
First of all, we remark that the flat solutions corresponding to \p = 3H 2 are 
reconstructed in the phase space irrespectively of the matter creation rate. 
In fact, taking the ratio of (|34]) to (|35|) we see that the parabola XP — 3H 2 



corresponds to the solutions for the flat case. Therefore, such a curve is a 
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separatrix for the general solutions even in presence of matter creation. This 
is just a reminder that (§) constrains open and closed spacetimes in such a 
way that they cannot evolve into one another. In general, there are at least 
two singular points H = p = in the phase plane: 



[H,p] = [0,0] , (36) 

[H,p] = [H d ,p d ] , (37) 

which correspond to Minkowski and de Sitter spacetimes. As will be shown 
ahead, the former plays a completely different role in this new context, while 
the later does not exist in the absence of matter creation. 

A phase diagram (or portrait) is a plot of the solutions to the dynamical 
system in the plane of the variables, with a flow corresponding to the arrow 
of time. As we know, this kind of graph can be compactified by changing 
the variables in order to bring infinity to a boundary in the phase portrait of 
the transformed dynamical system. A particularly useful new set of variables 



may be defined as the conformal mapping jl8 



H = H d ^— sin0 (38) 
1 — r 



XP = ^H 2 d [^—cos^ , (39) 

where H d is an arbitrary scale that will be set equal to the previously defined 
Hubble parameter of the de Sitter spacetime. The inverse transformation is 
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and 



easily computed as 



Ha 



(40) 



arctan ■ 



3H 2 

XP 



(41) 



It is readily seen that infinite values of p and H are brought onto the 
circle of radius r = 1. In terms of the new variables {r, 0}, we can recast the 
dynamical system (|54])-(PB]) as 



rHA sin (Acqs! 



37 



cos </> sin 2(p 



1 - 



"inHd sin < 



and 



(42) 



1 — r 



-Hh cos <t> cos ! 



1-| 



1-r) 



3nHd sin < 



(43) 



The phase portrait of usual FRW models (tp = 0; see, e.g., Ref. 0) is 
conformally transformed into the diagram of Fig.l, which have been shown 
for further comparison. The parabola XP — 3H 2 in the original phase space 
is mapped onto the radiuses at = ±7r/4, and the Big Bang (Big Crunch) 
singularities, respectively (p, H) = (oo, ±oo), are tamed to the points r = 1, 
cf) = ±7r/4. The shape of solutions is kept intact near the origin of the 
transformed graph. As expected, for these simple FRW models we are able 
to find exact expressions to the integral curves shown in Fig.l, since we know 
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these systems to be integrable (see, e.g., Ref. p9[). By taking the ratio of 
fl42|) to (]43| ) with -0 = we immediately solve the equation dr/d<p to find the 
solution 



r(0) = 1 - (1 - r ) 



,_2Z 

COS 0| I COS 20| 4 



2-37 

(44) 



where r Q and r are constrained to the interval [0,1]. Note that this formula is 
valid even for the unphysical region p < (cos < 0) which is also depicted 
for completeness. This will later prove helpful in order to analyze the nature 
of equilibrium points such as the Minkowski spacetime at the origin r = 
(p = H = 0), and the physical singularities at r = 1 and = ±7r/4 (Big 
Bang and Big Crunch respectively). 

At this point we should recall that although in flat spacetimes the energy 
density and Hubble parameter are indistinguishable, in curved spacetimes 
this is not so. In this way, the ansatz proposed in Refs.|| §], -0fl a t = 3(3H 2 , 
must be somewhat generalized to include the presence of curvature. In prin- 
ciple, if we mimic the bulk viscosity mechanism, a more reasonable phe- 
nomeno logical law for \1/ should be a power- law dependence \& = rjp u , where 
7] is a dimensional constant and < v < 1 (see, for instance, Ref.|^]). In 
this case, regardless of the curvature parameter, the Prigogine et al. ansatz 
corresponds to v = 1 and rj = f3\- However, this kind of phenomenological 
law does not include the case \l/ = 3j3nH {(3 constant), discussed in Ref.[[l^ . 
In this connection, we recall that the matter creation rate is a degree of free- 
dom introduced in the theory trough the balance equation for the particle 
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number density. In this way, we assume here a rather general expression 
which include both cases, namely: 

^ = 3^n(^)"V +1 . (45) 

where (3 is a constant. 

As noted above, a general feature of models with matter creation is the 
presence of a de Sitter solution. It is easy to see how these solutions come 
up: consider the curve C(p, H) corresponding to the equality ijj(p, H) = 3nH 
which, from Eqs.(|33[)-(|3SJ) always crosses the parabola \P = 3H 2 (or, equiv- 
alently, the straight line = n/4 of Fig.l). At that precise point [Hd, pd] 
(hereafter indicated as X), it is readily seen that both p and H are zero, 
meaning an exponentially expanding universe with constant density, i.e., an 
inflating universe. The question of whether models evolve towards it or out- 
wards from it is related to the nature of the singular point X. This is done by 
finding the Liapunov coefficients of the linearized dynamical system at that 
point ||16fl . For the matter creation rate of the form assumed above, after 
some simple algebra we find 

Ai = -2 

and 



A 2 = y(/i + ^) • 

This means that the fixed point X is either a saddle point, an attracting 
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Jordan node or a simple attracting node (attractor), if A2 > 0, A2 = 
or A2 < respectively. The absence of complex Liapunov coefficients for 
any 7 or (3 guarantees no closed orbits around X, which is obvious given 
the constraint (0) and the discussion after Eq.fl35|). The reader may verify 
for himself that A2 is negative if we consider ^ = j3p, since in this case 
v — — 1 and ^ = 2(7 — l)/7< 1. It thus follows that the inflating universe 
is an attractor for the class of models considered in Refs.|| |9|. The phase 
portrait is shown in Fig. 2. As explained, the point X is an attractor, and 
all expanding solutions converge there. Since H = (Ji. in Fig. 2) is a line 
of singular points, expanding and contracting models are never in contact. 
Such an effect is probably due to the creation pressure p c which diverges at 
that line (see Eq.([7|)). Actually, it is easy to see that all along line Ti the 
rate H (see Eq.(|3^)) goes to infinity. 

In this model, the expanding universe seems to have appeared either from 
a FRW-type singularity at B (r = 1, <fi — 7r/4), or at any point in the singular 
line Ti. In other words, all expanding spacetimes evolved from a singularity 
of some kind. Even Minkowski (denoted Ai), which is a nonsingular saddle 
point in FRW models (see fig.l), is a multiple equilibrium singular point here. 
Note also that closed universes show no turning points, either expanding or 
contracting forever. 

Contracting universes in this scheme can begin either from Minkowski or 
from the same singular line as the expanding ones, 7i. All solutions now 
evolve to a crunch, exactly as in the FRW case. In the region of negative 
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energy (left half of Fig. 2), contracting universes emerge either from the sin- 
gular line TL or from A4, and finish with infinite negative energy at (r = 1, 
= — 7r). Expanding universes with negative energy all start from Ti, and 
end up as flat empty universes at M.. In sum, if the universe described by 
this model started from a fluctuation with both H > and p > , it would 
evolve towards the inflating universe at X. On the other hand, if it started 
from a curvature fluctuation (H ^ 0, p — 0), it would certainly end up ei- 
ther in a crunch or in an empty Minkowski spacetime. The case of a matter 
fluctuation (H = 0, p ^ 0) is somewhat ill-defined, since p c diverges on the 
line H. 

5 Conclusion 

In this paper, by allowing both pressure and curvature, we have generalized 
the flat pre-inflationary stage with matter creation proposed by Prigogine 
and coworkers. All expanding solutions evolve to the flat de Sitter spacetime 
in the infinite cosmic time of their evolution. Physically, open and closed 
models are always singular, while the flat solutions split into two distinct 
subclasses which depend on the matter creation rate. The beginning of the 
universe in this case may be either a Minkowskian vacuum or a FRW type 
singularity. The avoidance of a singularity depends on the strength of the 
matter creation process, which has been translated here as a condition on the 
dimensionless S parameter. It thus appear that a certain degree of fine tuning 
is needed in order to escape from a physical singularity in these models. 
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Captions for Figures 



Fig.l 



Fig. 1 - Phase diagram of the standard FRW models mapped onto a 
disk. Flat universes are represented by the straight lines at <j) = ±7r/4. 
The origin M is an attracting node. Point B corresponds to the Big 
Bang singularity, while C is the Big Crunch at the end of the evolution 
of closed universes. 



Fig. 2 



Fig. 2 - Phase diagram of a model with ip = PxP- Again flat universes 
corresponds to straight lines at <p = ±7r/4. However, the point X is an 
attractor for all expanding solutions. H (corresponding to H = 0, or 
4> = 0) is a line of singular points and Minkowski (M) now is a Saddle 
point. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9606064v 1 



This figure "figl-2.png" is available in "png" format from: 
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